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Pu ABSTRACT 

We present deep broad-band imaging and long-slit spectroscopy of three compact, low-mass 
starburst galaxies at redshift z '^0.2-0.3, also referred to as Green Peas (GP). We measure physi- 
cal properties of the ionized gas and derive abundances for several species with high precision. We 
find that the three GPs display relatively low extinction, low oxygen abundances, and remark- 
ably high N/0 ratios We also report on the detection of clear signatures of Wolf-Rayet (WR) 
stars in these galaxies. We carry out a pilot spectral synthesis study using a combination of 
both population and evolutionary synthesis models. Their outputs are in qualitative agreement, 
strongly suggesting a formation history dominated by starbursts. In agreement with the presence 
of WR stars, these models show that these GPs currently undergo a major starburst producing 
between ~4% and '^20% of their stellar mass. However, as models imply, they are old galaxies 
having had formed most of their stellar mass several Gyr ago. The presence of old stars has been 
spectroscopically verified in one of the galaxies by the detection of MgiAA5167, 5173 absorption 
line. Additionally, we perform a surface photometry study based on HST data, that indicates 
that the three galaxies posses an exponential low-surface brightness envelope. If due to stellar 
emission, the latter is structurally compatible to the evolved hosts of luminous BCD/Hii galaxies, 
suggesting that GPs are identifiable with major episodes in the assembly history of local BCDs. 
These conclusions highlight the importance of these objects as laboratories for studying galaxy 
evolution at late cosmic epochs. 

Subject headings: galaxies: abundances — galaxies: dwarf — galaxies: evolution — galaxies: starburst 

1. INTRODUCTION vey (SDSS) Data Release 7 by volunteers in the 

"Galaxy Zoo" project (Lintott et al. 2008, 2011). 
"Green Peas" (GPs) are convenient laborato- Their nickname reflects their point-like appear- 

ries to study galaxy assembly at relatively low ^nce and "green" color on SDSS image over- 

redshifts (0.11 <z< 0.35). These compact, high- ^^ys. The latter is a consequence of very strong 

surface brightness systems were recently identi- [Olil]A5007 line emission, with equivalent widths 

fied on images from the Sloan Digital Sky Sur- _2000 A, enhancing the observed SDSS 

r fluxes at those redshifts. Such an extreme neb- 

iCONSOLIDER-GTC fellow; ^^a^. e^iiggion contribution has so far been doc- 
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umented m extremely metal-poor blue compact 
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dwarf galaxies (XBCDs) only (Tcrlcvich ct al. 
1991; Izotov et al. 1997; Papaderos et al. 1998, 
2008) in addition to a few ultra-compact star- 
bursting dwarfs in galaxy clusters (Reverts et al. 
2007). 

The GPs were first studied in detail by Carda- 

monc ct al. (2009) (hereafter C09), who showed 
that these galaxies reside in lower-density envi- 
ronments and are very rare (^2 galaxies deg^^ 
brighter than 20.5 mag). A subset of 80 GPs 
with decent signal-to-noise (S/N) SDSS spectra 
was spectroscopically characterized as purely star- 
burst systems. On average, these galaxies appear 
luminous in both optical {Mb ~ — 20 mag) and 
UV {Lpuv ~ SxlO^'^L©) wavelengths, and are 
characterized by high surface brightness and very 
compact appearance (typical sizes <5 kpc). Ac- 
cording to C09, the GPs are low-mass galaxies 
(stellar masses < 10^°"''Mq) with prodigious 
star formation rates (SFR up to 60 yr^^) and 
low intrinsic reddening (E(B-V)<0.25). In par- 
ticular, their specific star formation rates (sSFR 
in the range 10"'' to 10~^ yr"^) are among the 
highest inferred in the nearby Universe (cf e.g., 
Brinchmann et al. 2004; Salim et al. 2007), and 
well in the range of those of high-redshift galaxies 
(e.g., Bauer et al. 2005). 

Chemical abundances of the ionized gas have 
provided important additional clues about the na- 
ture of the GPs. Oxygen abundance determina- 
tions, based on the direct (Te) method led Amorfn 
et al. (2010) (hereafter AlO) to conclude that the 
GPs in the Cardamone's sample are genuinely 
metal-poor galaxies, spanning a range of values 
7.6 < 12+log(0/H) < 8.4 with an average of one 
fifth of the solar value in their gas-phase nietallic- 
ity in their gas-phase metallicity. Their results 
were recently confirmed by Izotov et al. (2011) 
(hereafter Izll), who pointed out that GPs are 
a subset of luminous compact galaxies showing 
chemical abundances similar to lower-luminosity 
Blue Compact Dwarfs (BCDs). 

Interestingly, the position of the GPs in the fun- 
damental relation between stellar mass and metal- 
licity (the mass-metallicity relation, MZR), and 
between B— band luminosity and metallicity (the 
luminosity-metallicity relation, LZR) appear sys- 
tematically offset (up to ~0.3 dex in the MZR) to 
lower abundances when compared with the bulk of 
local star-forming galaxies (SFGs) from the SDSS 



(AlO, Izll, see also Amorfn et al. 2011). Their 
location in the MZR, and also in the LZR, appear 
to form a distinct sequence, along with nearby 
XBCDs (I Zw 18 or SBS 0335-052, Guseva et al. 
2009) and some luminous BCDs at low (Bergvall 
& Ostlin 2002, Izll) and intermediate redshifts 
(e.g., Hoyos et al. 2005; Kakazu et al. 2007; Salzer 
et al. 2009) and most SFGs at high redshift (e.g., 
Pettini et al. 2001; Erb et al. 2006; Perez-Montero 
et al. 2009; van der Wei et al. 2011; Finkelstein et 
al. 2011). 

In this context, ionic abundance ratios between 
species with an assumed different stellar origin 
are important to probe the chemical evolution of 
the GPs. This is the case of the nitrogen-to- 
oxygen ratio (N/0), since the nitrogen and oxy- 
gen yields are driven by stars of different mass, 
therefore it gives relevant information about the 
SF rate and history of SFGs (e.g. MoUa et al. 
2006). In metal-poor SFGs (12-hlog(0/H)<8), ni- 
trogen production is expected to have mainly a 
primary origin, owing to massive stars. Then, in 
the O/H vs. N/0 diagram, they form a plateau 
at log(N/0) ~ -1.5 (Alloin et al. 1979; Camp- 
bell et al. 1986; Izotov & Thuan 1999; Pilyugin 
et al. 2003), with a relatively large vertical dis- 
persion (e.g., Garnett 1990; Pilyugin 1993; Henry 
et al. 2006; Nava et al. 2006; van Zee & Haynes 
2006; Perez-Montero & Contini 2009). For ex- 
ample, BCDs in the range 12-hlog(0/H)~7.6-8.2, 
generally show low N/0 ratios between -1.54 and 
-1.27 (Nava et al. 2006), with few exceptions only 
(see also Fig. 1 in Henry et al. (2006)). By con- 
trast, the nitrogen production in metal-rich SFGs 
(12-|-log(0/H)>8.2) has mainly a secondary ori- 
gin, powered by low- to intermediate-mass stars, 
which produce a positive correlation between N/0 
and O/H (e.g., Vila Costas & Edmunds 1993). In- 
triguingly enough, GPs show systematically larger 
N/0 ratios compared to m,ost SFGs a,t the same 
oxygen abundance (AlO, see also Amorfn et al. 
2011; Pilyugin et al. 2012), being in most cases 
located above the plateau. Their N/0 ratios, how- 
ever, seem to be normal for SFGs of similar stellar 
masses. 

The known global properties of the GPs sup- 
port the view that these galaxies go through a 
short and extreme phase in their evolution. How- 
ever, in order to confirm this and elaborate a co- 
herent evolutionary picture for GPs, important 
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pieces of the piizzle need to be siipplicd and inves- 
tigated. One example is the star formation history 
(SFH) of GPs, which still is not well-constrained, 
to a large extent because quantitative studies on 
the photometric structure and mass contribution 
of an underlying older stellar host are lacking. 
In order to establish a physical and evolutionary 
connection between GPs and nearby BCDs, it is 
crucial to verify that the structural properties of 
the host galaxy in both SFG classes are compati- 
ble. Whereas such a structural similarity has been 
demonstrated for luminous compact blue galaxies 
(LCBGs) at z - 1 by Noeske et al. (2006), no sur- 
face photometry studies for the lower- 2; GPs exist 
as yet. Another example is the interpretation of 
the relations between oxygen and nitrogen abun- 
dances, and with stellar mass and SFRs. They 
were discussed by AlO in terms of the balance be- 
tween inflows of metal-poor gas (e.g., Koppen & 
Hensler 2005) and the presence of enriched out- 
flows (e.g., van Zee et al. 1998). However, other 
scenarios usually invoked to explain large N/0 ra- 
tios in nearby metal-poor SFGs, such as possible 
pollution by Wolf-Rayet (WR) stars (e.g., van Zee 
et al. 1998; Brinchmann et al. 2008; Monreal-Ibero 
et al. 2010; Berg et al. 2011) were not addressed. 

The main reason why the above questions re- 
main open is that, until now, GPs have been 
mainly studied using SDSS data. This has al- 
lowed to infer relevant global properties for a large 
number of galaxies. The downside, however, is 
that SDSS spectrophotometric studies are limited 
by the often poor signal-to-noise (S/N) and sen- 
sitivity of the data. These limitations prevent, 
for example, studies of faint spectral features, like 
those due to young WR stars or of weak stel- 
lar absorption features owing to an old underly- 
ing component. Moreover, the moderate to poor 
S/N of SDSS spectra impacts accurate chemical 
studies requiring precise measurements of, e.g., 
the [Oiii]A4363A or [Nii]A6584A line fluxes (e.g., 
AlO) as well as stellar mass determinations based 
on a refined modeling of the spectral energy distri- 
bution (SED) of the stellar and gaseous continuum 
(Izll). 

In this study, we aim at going one step fur- 
ther in the understanding of the properties of GPs 
by using deep imaging and spectroscopy with the 
10.4m Gran Telescopio Canarias (GTC), in addi- 
tion to archival HST images. These high-quality 



data allow us to analyze in better detail the chem- 
ical and structural properties of a small sample of 
GPs. In particular, deep spectra with the GTC 
arc used both to improve on chemical abundance 
determinations and to study the SFH of GPs us- 
ing population and evolutionary spectral synthe- 
sis models. Faint features, e.g., absorption lines 
from old stellar populations or signatures of mas- 
sive and young WR stars, are examined. The lat- 
ter will provide strong constraints on the age of the 
young starburst and on the amount of metal pol- 
lution in the interstellar medium that can be ex- 
pected from them as well (e.g., Schaerer & Vacca 
1998; Guseva et al. 2000; Brinchmann et al. 2008; 
Perez-Montero et al. 2011). 

The paper is organized as follows: Section 2 de- 
scribes the target selection, the GTC-OSIRIS ob- 
servations and the data reduction, and the HST 
data. In Section 3 we show the methodology and 
present the results, which are then discussed in 
Section 4. Finally, conclusions are given in Sec- 
tion 5. We assume a standard cosmology with 
Hq = 70, Qa = 0.7, and = 0.3. 

2. DATA 

2.1. TEirget selection 

For the present study we have selected three 
GPs, SDSS J004054.32-K153409.6, SDSS J113303.79-h651341.3, 
and SDSS J232539.22-F004507.2 (for simplicity 
hereafter GP004054, GP113303, and GP232539 
respectively) from C09. The main properties of 
the galaxies as derived from the literature are 
summarized in Table 1. The three galaxies are 
at very similar rcdshifts z ^ 0.24 — 0.28 and they 
are located in relatively isolated regions, with no 
nearby bright companions. These GPs have a high 
UV surface brightness and luminosity {^\0^^-^Lq 
), translating into large SFRs (> lOM© yr~^) 
per unit area. For this reason, they were also in- 
cluded in the sample of ^30 super compact, UV 
luminous galaxies by Hoopes et al. (2007), also 
known in subsequent studies as local "Lyman- 
break Analogs" (LBAs Basu-Zych et al. 2007; 
Overzier et al. 2008, 2009, 2010; Gongalves et al. 
2010). They were, therefore, better studied than 
other GPs. For example, these galaxies are the 
three out of four GPs for which spatially resolved 
HST imaging is publicly available up to date. 
Prom HST WFPC2 and ACS imaging (Overzier 
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et al. 2009) wc know that they share similar mor- 
phologies, with UV and optical light dominated by 
few very luminous star-forming clumps superim- 
posed on a compact stellar host (optical half-light 
radius of about 1 kpc). 

2.2. Data set 

Deep broad-band imaging and long-slit spec- 
troscopy for the three target galaxies were car- 
ried out using the OSIRIS instrument, mounted 
on the 10.4 m GTC at the Observatory Roque de 
los Muchachos (La Palma, Spain). 

OSIRIS^ (Ccpa et al. 2000) is an imager and 
spectrograph for the optical wavelength range 
(from ~3650 to lOOOOA), located in the Nasmyth- 
B focus of GTC. It consists of two 2048x4096 
Marconi CCD42-82 with a 9.2 arcscc gap between 
them. The unvignetted instrument field of view is 
7.8x7.8 arcmin with a pixel scale of 0.125 arcsec. 
Both imaging and spectroscopic observations were 
obtained in service mode by the GTC staff during 
the first semester of 2010. The log of observations 
is summarized in Table 2. 

2.2.1. GTC and HST Imaging 

Broad-band imaging for the three targets were 
obtained using the SDSS z' filter (centered at 
9695 A) in the standard mode (2x2 binning), given 
a pixel scale of 0.25 arcsec. In order to avoid 
substantial contamination from strong emission 
lines and reach faint surface brightness levels, we 
chose the z' filter despite its efficiency is lower 
than the i' filter. For each of the three galaxies 
our broad-band images were taken under differ- 
ent sky transparency conditions. Unfortunately, 
only GP113303 was observed in a dark night un- 
der photometric conditions. The seeing was in all 
the cases below 1.2 arcsec. Total exposure times of 
2250 sec were achieved taking several scries of five 
short exposures (90 sec) in a cross-shape dithering 
pattern with offsets of 10 arcsec. Series of bias, 
twilight sky flats, and several spcctrophotomctric 
standards were observed during the same nights 
(see 2.2.3). 

The main goal of deep imagery with the GTC 
was the investigation of the close environment of 



Detailed information on GTC and OSIRIS can be found in 
http://www.gtc.iac.es 



OUT sample GPs e.g., extended tidal low-surface 
brightness features that could have gone unde- 
tected on the shallower SDSS and archival HST 
images. In Figure 1 wc present deep z'— band im- 
ages for the three GPs. The right gray scale bar 
and contours in these images show the z'— band 
surface brightness of the galaxies. Only GP232539 
appears to show a somewhat extended (size > 
FWHM) LSB component. Possible small compan- 
ions, barely seen on SDSS images, are projected 
few arc seconds to the galaxies. 

In order to study the morphology and the struc- 
tural properties of our sample galaxies, we ad- 
ditionally included in our analysis archival HST 
WFPC2 images in the filter F606W (HST pro- 
posal ID 11107, P.I: T. Heckman), which are also 
presented in Figure 1 as insets. These images re- 
veal the complex morphology of the inner, high 
surface brightness regions of these GPs, and will 
be used for discussion in Section 4.2. 

2.2.2. GTC spectroscopy 

Long-slit spectroscopy was carricid out in the 
standard mode (2x2 binning) with the highest res- 
olution mode available at these dates, R~1018 (at 
5510A) and R-1122 (at 7510A). Thus, we used 
the RIOOOB and RIOOOR grisms and a slit width 
of 0.8 arcsec, projecting onto a full width at half 
maximum (FWHM) of about 3 pixels. This setup 
yields wavelength coverages in the blue ~3630- 
7500A and in the red ~5100-10250A, with typi- 
cal dispersion values of 2.1 and 2. 6 A pixel" at 
55 10 A and 75 10 A respectively. The spectra were 
taken along the parallactic angle. For each galaxy, 
blue and red spectra were taken during different 
(dark) nights. Seeing conditions varied between 1 
and 1.4 arcsec, while atmospheric conditions were 
also not uniform, being two of the nights spec- 
troscopic, three of them clear, and the remain- 
ing night slightly cloudy. Some Saharian dust was 
present in the atmosphere, affecting those obser- 
vations at lower air masses significantly, reducing 
their final S/N, and doing more difficult the sky- 
background subtraction. This is more evident in 
the red part of the spectra, were large sky resid- 
uals are present. Noisiest regions in the blue and 
red ends of the observed spectral range were not 
considered in the subsequent analysis. A series of 
bias, twilight sky and dome ffats, two calibration 
lamps, as well as one or two spectrophotometric 
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standards were observed during the same nights. 

2.2.3. Data reduction 

The data were fully reduced and calibrated us- 
ing IRAF^ routines. This includes the usual pro- 
cedures for bias and overscan subtraction, flat- 
fielding corrections, cosmic ray removal, and co- 
addition. For the broad-band images, additional 
large scale illumination and fringing patterns were 
removed effectively after subtraction of a "mas- 
ter sky" frame. This was obtained combining the 
5x5 dithered science frames (after standard cor- 
rections and masking-out all saturated stars in 
the field) of each object. Finally we got deep 
(/Uz~26-27 mag arcsec^^) images after a good sky- 
background subtraction. For the spectra, wave- 
length calibration was done using HgAr+Xe+Ne 
lamp arcs. The accuracy (<0.lA) was checked 
a posteriori using sky emission lines. The spec- 
tra was corrected for atmospheric extinction and 
then flux calibration was performed using sev- 
eral spectrophotometric standards (Grw70-|-8247, 
Ross640, L1363-3, G157-34) along the observing 
blocks, taken at the same night for each grating. 

3. GTC SPECTROSCOPY: ANALYSIS 
AND RESULTS 

In Figure 2 we present the OSIRIS-GTC spec- 
tra for the three GPs. The spectra are dominated 
by intense narrow nebular emission lines on top a 
faint blue stellar continuum lacking stellar Balmer 
absorption features. The high S/N of spectra al- 
lowed for a detection of several faint emission lines 
(e.g. the temperature-sensitive [Olll]4363 line) 
that are displayed magnified in the lower panel. 

3.1. Emission-line intensities and redden- 
ing correction 

We measured emission line intensities on the 

reduced and calibrated spectra of the three galax- 
ies using the task splot of the package iraf. To 
measure the flux of a given line, we integrated the 
flux between two points given by the position of 
a local continuum placed by eye. The statistical 



^iraf: the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observato- 
ries, which is operated by the Association of Universities 
for Research in astronomy, Inc. (AURA) under coopera- 
tive agreement with the National Science Pundation 



errors associated with the measured emission line 
fiuxes were calculated using the following expres- 
sion (Perez-Montero & Dfaz 2003): 

(7j = (Tcy ^ + 

where di is the error flux of the measured emis- 
sion line, fJc is the rms error derived for the lo- 
cal continuum, A'^ is the number of pixels taken 
for the measurement of emission line flux, Wi is 
the absolute value of the emission line equivalent 
width, and A is the wavelength dispersion. Al- 
though there is absorption of the Balmer emission 
lines caused by underlying stellar populations in 
the objects (Dfaz 1988), we checked in the resid- 
uals to the STARLIGHT fitting (see § 3.3) that this 
effect is negligible compared to the reported er- 
rors. The emission line fluxes, F(A), relative to 
F(H/3) = 1000, are listed with their corresponding 
errors in Table 3. 

Each emission line flux is affected by the pres- 
ence of interstellar dust which absorbs it according 
to the law: 

I{HI3) F{Hp) 

where 1(A) and F(A) are the corrected and mea- 
sured emission line; fluxes, respectively, c(H^) is 
the constant of reddening, and f(A) is the ex- 
tinction law, which for this work, we took from 
Cardelli, Clayton, & Mathis (1989), and whose 
values are also listed in Table 3 for the correspond- 
ing listed emission lines. c(H/3) was calculated 
for each object as the error- weighted least square 
flt to the relation between the extinction law and 
the quotient between the observed-to-theoretical 
Balmer and H/3 emission line fluxes before and af- 
ter underlying continuum subtraction. We found 
that the differences between the corrected and un- 
corrected c(H/3) values are always smaller than 
the uncertanties quoted in Table 3. The Balmer 
emission-lines with enough S/N reach up to H17 
for GP004054, Hll for GP1133G3, and H15 for 
GP232539. The theoretical emission-line ratios 
were derived using Storey & Hummer (1995) val- 
ues for the appropriate electron density and tem- 
perature for each one of the three galaxies. Tak- 
ing the three objects coordinates into account to 
calculate the corresponding Galactic extinction, 
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which arc Cg(H/?) 0.09 for GP004054, 0.02 for 
GP113303, and 0.05 for GP232539, we conclude 
that the dominant contribution to reddening is due 
to intrinsic extinction. The extinction-corrected 
emission hne fluxes relative to I(H/3) = 1000, along 
with their corresponding errors are listed in Table 
3. In this tabic we also list the respective con- 
stants of reddening, the extinction corrected H/3 
flux and the H/3 equivalent widths. 

3.2. Electron densities and temperatures 

All physical conditions, including electron den- 
sity and temperature, were calculated using the 
task temden of the package iraf, taking for each 
species the same atomic coefficients as in Hagele 
et al. (2006). All of them, as described in this 
subsection, are listed in Table 4. 

Electron density has been estimated by taking 
the emission-line ratio of [Sii] 6716, 6731 A. The 

error propagation from the emission-line fliixes 
does not allow to give a precise estimate of the 
density in any of the three galaxies, but it gives in 
all cases an upper limit which is much lower than 
the critical density for collisional deexcitation. 

Electron temperature of [Oiii] was derived with 
high precision in the three galaxies by taking the 
emission line ratio between the sum of 4959, 5007 
A and 4363 A. This gives temperatures for the 
three objects which are typical for BCDs and 
Hii galaxies (Campbell et al. 1986; Masegosa et 
al. 1994; Izotov et al. 2006; Kehrig et al. 2006; 
Hagele et al. 2006, 2008), ranging from 13400 K 
for GP004054 and 14600 K for GP232539. 

Other electron temperatures than [Om] were 
not derived directly in any of the three objects. Al- 
though [Oil] 7319, 7330 AA and [Nil] 5755 A were 
detected in the spectra of the galaxies, they did 
not lead to a precise derivation of the correspond- 
ing electron temperatures. Therefore, in order to 
calculate the corresponding ionic abundances, we 
derived these temperatures from te([Oiii]) and us- 
ing the expressions obtained from photoionization 
models described by Perez-Montero & Diaz (2003) 
and Perez-Montero & Contini (2009) for te([Oll] 
and te([Nii]), respectively. Regarding te([Sill]), 
as the nebular lines at 9069, 9532 A have not 
been detected, we derived this temperature again 
from te([Oiii]), but using the empirical relation 
obtained by Hagele et al. (2006). These derived 



temperatures are also listed in Table 4. 

3.3. Ionic and total chemical abundances 

Ionic abundances were calculated for the visible 
chemical species in our optical spectra by using 
the task ionic of the IRAF package and taking 
the same atomic coefficients described in Hagele 
et al. (2006) and the corresponding ionization cor- 
rection factors (IGF). The ionic and total abun- 
dances, along with these ICFs are listed in Table 4. 

Helium abundances^ have been obtained as- 
suming that 

He _ He+ + He^+ 
IT ~ W 

He"*" abundances (noted as y"*" in Table 4) were 
calculated as the error- weighted mean of the abun- 
dances derived using the Hei emission lines at 
4471, 5876, 6678 and 7065 A, with te([Oiil]) and 
taking the expressions described by Olive & Skill- 
man (2004). Only in the case of GP113303 we 
measured the Hen nebular emission line, but the 
corresponding y^+ abundance is negligible as com- 
pared to y+. The derived He abundances (denoted 
as y in Table 4) range from 0.087 to 0.089. These 
values fall between the pre-galactic He/H ratio of 
0.08 (Luridiana et al. 2003) and the solar ratio 
0.10 (Lodders 2003), and are "normal" for low- 
metallicity star-forming galaxies (e.g., Kobulnicky 
& Skillman 1996; Hagele et al. 2008). 

Oxygen abundances were calculated assuming 
that 

O _ 0+ + 

with 0+ abundances derived from the [On] 3727A 
emission line and taking te([Oii]). In the case 
of 0^+, we toke [Oiii] 4959, 5007 A emission 
line intensities with te([Oiii]). The oxygen total 
abundances, ranging from 12+log(0/H) — 7.81 
for GP232539 to 7.98 for GP004054, are ~ one 
fifth solar [12-hlog(O/H)0 = 8.69, Asplund et al. 
(2009)]. 

Nitrogen abundances were calculated by assum- 
ing: 



'No correction by neutral Helium nor by collisional excita- 
tion have been considered 
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N 
H 



ICF{N+) 



N+ 
11+ 



and deriving N+ abundances from the [Nii] 6548, 
6584 AA emission-line fluxes and te([Nii]). The 
ICF was calculated using the expression described 
in Perez- Montcro & Contini (2009). The derived 
N/0 ratios, ranging from log(N/0) = -0.99 for 
GP232539 to -1.06 for GP004054, arc well above 
the typical values corresponding to star-forming 
galaxies of the similar metallicity, i.e., the N/0 
plateau seen at log(N/0)~ —1.5, (e.g., Izotov & 
Thuan 1999). 

Neon abundances were calculated by assuming: 



^ « ICF{Ne^+) 
H 



H+ 



calculating Ne^+ ionic abundances from [Neiii] 
3868 A with te([Olll]), and using the expression 
for the corresponding ICF described by Perez- 
Montcro ct al. (2007). The derived Ne/0 ra- 
tios arc all slightly higher than the solar value 
[log(Ne/O)0 = -0.76, Asplund et al. (2009)]. 

Argon chemical abundances were calculated in 
GP004054, for which both [Arm] 7136 A and 
[Ariv] 4740 A were measured, and in GP232539, 
with the [Ann] 7136 A emission line. Then, since 
we have both ionic argon abundance Ar^+ and 
Ar3+ in GP004054, and only the first one, Ar^^, in 
GP232539, different ICFs were considered accord- 
ing to the expressions proposed by Perez-Montero 
et al. (2007). As in the case of neon, both argon- 
to-oxygen ratios result higher than the solar value 
[log(Ar/O)0 = -2.29, Asplund et al. (2009)]. 

Sulphur abundance was calculated by assuming 
the following: 



^ = ICF{S+ + 52+) • 
H 



S+ + 
H+ 



with the S+ derived from the [Sii] 6717, 673lAA 
and assuming that te([Sii]) ~ te([Oii]). In the 
case of , since our available spectral range has 
not allowed a measurement of the [Siii] emission 
lines at 9069, 9532 AA, we derived their ionic 
abundances with [Slll] 6312 A, and with te([Slll]). 
We considered the ICF for (S+-t-S2+) obtained by 
Perez-Montero et al. (2006). The derived sulfur- 
to-oxygen ratios range from -1.22 for GP113303 



to -1.55 for GP232539. The high quoted errors 
for S/0 make the derived values to be consistent 
with the solar value [log(S/O)0 = —1.57, Asplund 
et al. (2009)]. 

Finally, iron abundance was calculated from 
Fe^"*" with the emission line relative intensity of 
[Feiii] 4658 A and the electron temperature of 
[Oiii]. We have also used the ICF proposed by 
Rodriguez & Rubin (2004). All Fe/0 are simi- 
lar in the three galaxies, ranging from -1.49 for 
GP004054 to -1.69 for GP232539. 

The derived Ne/0, S/0, Ar/0, and Fe/0 val- 
ues (see Table 4) are, within uncertanties (which 
are especially high for S/0 and Fc/0), consistent 
within the three galaxies. Overall, their values 
can be considered as "normal" when compared to 
the mean values found for nearby BCDs and Hll 
galaxies in the literature (e.g., Izotov & Thuan 
1999; Hagele et al. 2006, 2008). 

3.4. Spectral fitting 

In order to gain insights into the SFH of the 
GPs under study, we interpreted their integrated 
spectra by means of spectral synthesis models. 
To check the consistency of our results, we used 
both the population synthesis code STARLIGHT'* 
(Cid Fernandes et al. 2004, 2005; Mateus et al. 
2006) and a two-component evolutionary synthesis 
code that is based on PEGASE 2.0 (Fioc & Rocca- 
Volmerange 1997). 

We used starlight to synthesize the observed 
stellar continuum of the galaxies as due to the 
superposition of single-age stellar populations 
(SSPs) of different ages and metallicities. We 
used the SSP library provided by the popstar^ 
synthesis models (hereafter run 1, MoUa et al. 
2009; Martfn-Manjon et al. 2010, Garcia- Vargas 
et al. in prep.). These models follow the evolution 
of SSPs from very young (0.1 Myr) to very old 
(15.8 Gyr) ages by combining Padova '94 stellar 
evolution models with the most recent physics for 
stellar atmospheres and nebular continuum emis- 
sion. POPSTAR models are therefore particularly 



*The STARLIGHT project is supported by the Brazil- 
ian agencies CNPq, CAPES and FAPESP and 
by the France-Brazil C APES/Cofecub programme, 
http:/ /www. starlight. ufsc.br. 

® POPSTAR models arc publicly available from 
http;//www.fr actal-es.com/PopStar/SEDmod. html 
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well suited for the modeling of star-forming galax- 
ies (Molla et al. 2009) . Note that the low spectral 
resolution (20 A) of the currently available pop- 
star SSPs docs not permit a perfect match to 
Balmer stellar absorption features. This, how- 
ever, is not expected to have a notable impact on 
the derived SFHs. 

In low-metallicity starburst galaxies, a signif- 
icant contribution from nebular continuum emis- 
sion is expected to be superimposed upon the stel- 
lar SED. The gaseous emission not only enhances 
the luminosity of the galaxy but it additionally 
tends to make its spectral continuum redder than 
what is expected from a stellar SED (e.g., Krueger 
et al. 1995; Papaderos et al. 1998, Izll). There- 
fore, as recently pointed out by Izll, if nebular 
continuum emission is not taken into account in 
SED fitting for starburst galaxies the estimated 
stellar masses can be severely overestimated. In 
this respect, the new popstar SSPs including 
nebular continuum and line emission offer an im- 
portant advantage towards a realistic SED model- 
ing of GPs. 

In modeling the stellar SED we used SSPs 
with three different metallicities, 0.008, 0.004, 
and 0.0004 (i.e., -1/2, 1/5, and 1/45 solar) for 
a Kroupa initial mass function (IMF) between 
0.15 and 100 Mq, and all available SSPs be- 
tween logt = 5.0-10.2 yr. For the SSPs with the 
youngest ages, i.e., <20 Myr, we only used mod- 
els with Z = 0.004, which is the closest value to 
the gas-phase metallicity that we determined (see 
Table 3) whereas for the older stellar component 
all three available metallicities were used. Prior 
to STARLIGHT models, the flux calibrated spec- 
tra were de-redshifted and resampled to 1 A/pixel. 
Spectral regions with strong emission lines or sky- 
subtraction residuals were masked out from fit- 
ting. Models were applied on the wavelength range 
3500-6700A (rest frame) to exclude noisy spectral 
regions. The S/N ratio of the spectra in the rela- 
tively featureless window between 4220 and 4280A 
varies from 15 to 45 for the three galaxies. 

3.4- 1- run 1: POPSTAR models 

The STARLIGHT fits based on the popstar 
SSPs are shown overlaid with the observed spec- 
tra in the upper panels of Figures 3-5. It can 
be seen that in all cases the synthetic stellar and 
gaseous SED provides a good match to the obser- 



vations. The right upper diagram shows the lu- 
minosity contribution (%) of the SSPs evaluated 
by STARLIGHT to the normalization wavelength of 
4170 A. Vertical thin-gray lines depict the ages 
available in the POPSTAR library. The stellar mass 
fraction corresponding to each SSP is plotted on 
the lower panel. 

The main output from the STARLIGHT mod- 
els for POPSTAR SSPs is summarized in Table 5. 
Columns 2-4 list the predicted relative contri- 
bution of the nebular emission to the total con- 
tinuum emission at the [On], II/3 and Ha wave- 
lenghts. The obtained fraction of young (<100 
Myr) stars M^^young with respect to the existing 
stellar mass of ~20%, suggests that GP004054 cur- 
rently experiences a major episode in its assem- 
bly history. The last two columns of Table 5 list 
the existing total stellar mass total and the re- 
duced X^. It is interesting to point out that the 
stellar mass M*. 

total estimated from the fits, '^2— 
3x10^ Mq, places GPs in the range of luminous 
BCDs (e.g. Bergvall & Ostlin 2002; Guzman 
et al. 2003), even when aperture corrections (cf. 
Sect. 3.4.4) are taken into account. It can be seen 
that the STARLIGHT solutions delineate two main 
SSP groups, a younger one dominating the optical 
light, and an older stellar component (~10 Gyr), 
contributing most (>80%) of the stellar mass. Sig- 
natures for an intermediate-age stellar population 
were found for GP232539 only, where SSPs at 
~0.25 and ~2.5 Gyr appear to make a substan- 
tial (> 20%) contribution to the stellar mass. 

Whereas the strong starburst activity in the 
GPs under study is obvious, already from their 
high emission-line EWs, the inferred mass fraction 
young of stars recently produced is to be consid- 
ered with some caution. Spectral synthesis models 
are known to be plagued by substantial degenera- 
cies, in particular for SFGs (see e.g. Guseva et 
al. 2001, for a detailed discussion) making quan- 
titative statements on the reality and relative im- 
portance of individual features in the derived SFH 
difficult. Additionally, it should be called into at- 
tention that here we present a pilot attempt of us- 
ing POPSTAR SSPs in conjunction with starlight 
and no rigorous tests of possible numerical effects 
in this context have been made so far. 

We therefore include below two further model- 
ing attempts in order to gauge the variation of 
the derived Mj,_young for the GPs under study. 
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In the first one (hereafter run 2), we employed 
STARLIGHT using, however, SSPs based on the 
stellar models by Bruzual & Chariot (2003). One 
of the main differences between popstar and 
Bruzual & Chariot models is the inclusion of neb- 
ular continuum emission in the former (a detailed 
comparison of popstar and other models can be 
found in Molla et al. (2009)). Additionally, the 
Bruzual & Chariot SSPs have a higher spectral 
resolution and assume a Salpeter IMF between 0.1 
and 100 Mq. 

3.4.2. run 2: B&C models 

For RUN 2 (lower panels of Figures 3-5) we im- 
posed roughly the same constraints as for RUN 1, 
using the same masks for emission lines (shaded 
areas) and normalization wavelength, and an SSP 
library with nearly the same age and metallic- 
ity coverage. The main results as derived from 
RUN 2 are also included in Table 5. Comparing 
them with those from RUN 1, we found a satis- 
factory agreement in the overall SSP age distribu- 
tion, with two prominent peaks at young and old 
ages for CP 11 3303 and signatures of an additional 
intermediate-age stellar population in GP004054 
and GP232539. In view of the differences in the 
SSP libraries and the details of the SFH obtained, 
some differences in, e.g., total stellar mass and 
young stellar mass fraction obtained from RUN 1 
and 2 are not surprising. With the exception of 
GP004054, RUN 1 yields a significantly lower total 
stellar mass and a higher mass fraction from young 
stellar populations by a factor of about 2 (~0.3 
dex) compared with those obtained from run 2 
(see Table 5). 

3.4.3. run 3: evolutionary synthesis models 

As a second consistency check (run 3), we 
applied a modified version of the evolution- 
ary synthesis code PEGASE 2.0 (Fioc & Rocca- 
Volmerange 1997). In this approach, each spec- 
trum was modeled as due to the superposition of 
the SED from an old and a young stellar popu- 
lation approximating, respectively, the underlying 
host galaxy and the starburst component. The 
host was modeled by an exponentially decreasing 
SFR since 13 Gyr with an e-folding time of 3 Gyr. 
Note that the assumed SFH for the host implies 
for GP004050 and GP113303 some contribution 
from intermediate-age to young (0.1-1 Gyr) stars 



which is not apparent from the starlight fits 
(run1&:2). As for the young stellar component, 
we assumed an instantaneous burst whose age was 
allowed to vary between and 50 Myr. The SEDs 
were computed assuming a fixed stellar metallicity 
of Z0/5 and a Salpeter IMF (0.1-100 Mq), and 
include nebular continuum and line emission. 

The best-fitting solutions were constrained by 
varying the burst age, extinction C(H/3) and bust 
parameter hpar (mass fraction of stars formed in 
the burst with respect to the mass of the stars ever 
formed) that reproduce best the observed SED 
continuum and Ha and H/3 EWs. We note that 
the concept used here, i.e. the exploration of the 
SFH and C(H/3) that self-consistently account for 
Balmer line EWs, in addition to stellar SEDs and 
colors of SFGs was originally used for specific tasks 
in Izotov et al. (1997) and Papaderos et al. (1998) 
and further developed in Guseva et al. (2001) and 
Guseva et al. (2007). 

In the lower panels of Figures 3-5 we show the 
best-fitting SEDs for the young and the old stellar 
component as blue and red shaded areas, respec- 
tively. It can be seen that the superposition of a 
burst on an evolved galaxy host (black curve) can 
reproduce quite well the spectra of GP004054 and 
GP232539, whereas some systematic deviations 
are present in the case of GP113303. This might 
be attributed to the absence of an intermediate- 
age population in the latter system, as revealed by 
the STARLIGHT models in RUN 1&2. The relevant 
model output is summarized in Table 6. Inter- 
estingly, in all cases the best-fitting models imply 
an C(H/3) very close to the measured value for the 
nebular component. They also yield a good match 
to the observed Ha and H/3 EWs. 

Tables 5 & 6 show the significant variation of 
certain parameters (e.g. M*, young) over run 1 
through RUN 3. While they may not be surpris- 
ing in view of the different models, SSP libraries 
and fitting constraints used, they reflect the in- 
herent uncertainties in the reconstruction of the 
SFH of individual starburst galaxies using state- 
of-art spectral synthesis models. On the other 
hand, it is worth pointing out that all models, es- 
pecially those including nebular emission (run 1 
and RUN 3), consistently suggest that our sample 
GPs presently undergo a significant evolutionary 
stage in which they rapidly form between a few % 
and up to -^20% of their stellar mass. This high- 
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lights the importance of these objects as labora- 
tories for studying galaxy evolution at late cosmic 
epochs. 

3.4-4- Aperture corrections 

We did not applied any aperture correction be- 
fore fitting our spectra. However, our flux cali- 
bration was checked using the SDSS spectra and 
SDSS broad-band magnitudes. We calculated 
the fraction of total optical light covered by the 
OSIRIS long-slit, deriving a slit coverage a factor 
of 1.5-1.8 smaller than those of the SDSS fiber. 
According to the fraction of light inside the SDSS 
fiber calculated in AlO, we estimated a slit cov- 
erage of about 45-55% of the visible light in the 
SDSS images. Therefore, we must keep in mind 
that all global properties that can be derived from 
the spectral analysis (e.g., stellar mass) may be 
lower than their true values. 

3.5. Wolf-Rayet features 

All massive (M > 25 Mq for Zq), very lu- 
minous (10"'' — IO^Lq) O stars pass through the 
WR phase 2—5 Myr after their birth, spending less 
than 5 X 10^ yr in this phase (e.g. Schaerer & Vacca 
1998, Meynet & Maeder 2005). The presence 
of WR stars, as often observed in extragalactic 
Hll regions (Gonzalez-Delgado et al. 1994; Garcia- 
Vargas et al. 1997), and in both integrated (e.g., 
Allen et al. 1976; Kunth & Sargent 1981; Conti 
1991; Masegosa et al. 1991; Vacca & Conti 1992; 
Schaerer et al. 1999; Guseva et al. 2000; Brinch- 
mann et al. 2008) and spatially resolved (Kehrig et 
al. 2008; Monreal-Ibero et al. 2010; Lopez-Sanchez 
et al. 2011) spectra of starburst galaxies, is char- 
acterized by the detection of two broad emission 
features in their optical spectra. One of these fea- 
tures is the blue WR bump, which is a blend of 
the Nlll A4640, Clll/CiV A4650, and the broad 
Hell A4686 emission line. This feature is mainly 
due to WN stars. In contrast, WC stars are the 
main responsible for the red WR bump, which is 
a blend of the Cm A5698 and Civ A5808 broad 
emission lines. The red WR bump is weaker than 
the blue one and is rarely observed at low metal- 
licities (Crowther 2007). WR stars are present in 
galaxies with recent starbursts (e.g., BCDs, Gu- 
seva et al. 2000; Perez-Montero et al. 2010). In 
CPs, however, no WR features had been detected 
so far. 



Figure 6 shows a zoom of the OSIRIS spectra 
in the 4200A-5200A (rest-frame) range, where 
the blue WR bump is detected in the three 
CPs. These WR features are also identified af- 
ter the subtraction of the continuum fitting by 
STARLIGHT. In contrast, none of the galaxies has 
a red bump detection. 

Given the low spectral resolution of our data, 
we measure the WR bumps by using the IRAF task 
splot. First, we fitted and substracted the adja- 
cent continuum. Then we fitted gaussians to the 
nebular emission lines in the region between 4600 
and 4750A, and finally we integrated the; remain- 
ing emission above zero to have an estimate of the 
blue WR bump fluxes. The main properties de- 
rived from these measurements are summarized in 
Table 6. For example, from the luminosity of the 
bump, and adopting a theoretical WR luminosity 
from STARBURST99 models (Leitherer et al. 1999) 
at the same metallicity, we inferred the number of 
WR stars expected for each galaxy. 

4. DISCUSSION 

4.1. Stellar mass 

In Table 5 we list stellar masses contained in the 
spectroscopic aperture as derived from the spec- 
tral fltting, and using the adopted distances to 
the galaxies (see Table 1), and mass-to-light ratios 
obtained from the mixture of popstar synthetic 
stellar populations fltted to each galaxy spectrum. 
We also give in Table 5 an estimate of the percent- 
age of the stellar mass given by stellar populations 
younger than 20 Myr, and the fractional contribu- 
tion — for three different spectral ranges — of the 
nebular continuum. 

Our stellar mass estimates from run 1 for 
GP004054 (M^ =1.7x1O9M0), GP113303 (Af* =2.8x109Mq), 
and GP232539 (M,, =2.4x1O3M0), confirm that 
these galaxies are dwarfish systems having had 
formed between a few % and up to ^-^20% of their 
present stellar mass in a recent (<100 Myr) ma- 
jor star formation episode. Note that the stellar 
masses derived here are not far from those in- 
ferred for the same objects in previous studies (see 
Table 1), despite significantly different modeling 
approaches used. 

Discrepancies are likely due to the inclusion of 
the contribution of the nebular continuum emis- 
sion (run 1) in this study. For the three analyzed 
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objects, we find this contribution to bo significant 
and to increase to several % towards the red spec- 
tral range (see Table 5). The impact of the red 
nebular continuum on SED fitting is already sug- 
gested by the fact that starlight models based 
on purely stellar SSPs (run 2) yield a by a fac- 
tor ~2 lower M^^^young, consequently a larger mass 
fraction of old stars. Precisely this effect was dis- 
cussed by Izll and invoked in order to explain the 
differences between their Mi, values and those de- 
rived by C09 for the whole sample of GPs. Note, 
however, that we do not find any clear correlation 
between the offset of the derived stellar masses 
and the relative flux contribution of the nebular 
continuum in our three studied objects. 

For GP004054, C09 and Izll did not published 
values. However, we compared our estimates 
with those given by Overzier et al. (2009). The lat- 
ter were obtained from the SDSS/DR7 stellar mass 
catalog ^, and derived using SDSS photometry and 
a large grid of models constructed on the basis of 
the Bruzual & Chariot (2003) library. Interest- 
ingly, stellar masses calculated here for GP004054, 
and also for GPl 13303 and GP232539 (also listed 
by Overzier ct al. (2009)) are within ^0.3 dex in 
agreement with our values. It is worth noting that 
values given by Overzier et al. (2009) did not take 
into account any contribution from nebular con- 
tinuum, which may be possibly masked or com- 
pensated by other uncertainties. 

4.2. The star formation history of the GPs 

Several constraints on the SFH of the studied 
GPs can be inferred from our results. As consis- 
tently indicated by spectral synthesis models, none 
of them presently forms its first stellar population. 
Quite to the contrary, an evolved stellar compo- 
nent with an age between a few 10* yr and several 
Gyr is present in all three GPs and provides at 
least ~80% of their stellar mass. In particular, di- 
rect evidence for an old stellar component is found 
for the system GPl 13303, where we detect a faint 
{EW ~ 1.5A) but clear MgiAA5167, 5173 broad 
absorption line. This feature is due to the presence 
of old stars and, as shown in Fig 7, it is clearly vis- 
ible by eye in the spectra. At least for GP113303, 
the detection of this feature gives us the confir- 



' Available at http://www.mpa- 
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mation of the presence of old stars. Whereas for 
GP004054 and GP232539 the spectral fitting led 
to similar results than for GP113303, absorption 
features in their spectra are probably too faint to 
be detected. 

Both population and evolutionary synthesis 
models including the nebular component (run 
1&3) imply, on the other hand, that the GPs 
studied here are experiencing a significant stage 
in their buildup, as they form ^4-20% of their 
total mass in an intense intermittent or prolonged 
starburst episode which, as suggested by the pres- 
ence of WR features is very recent or still ongoing. 

The strong starburst (SFR~ 7 - IAMq yr'^ 
see Table 1) taking place in these galaxies and 
their low stellar masses (~ 1.5-2.7xlO^M0 yr~^) 
imply large specific star formation rates, sSFR~ 
2 — 9xl0~® yr~^. In agreement with C09 and 
Izll, these values are typical for GPs but they are 
unusually high in the local Universe (Brinchmann 
et al. 2004), being more comparable to those ob- 
served in high redshift starburst galaxies. Since 
these sSFRs imply relatively short mass doubling 
times (= 1/sSFR) <500 Myr), it is likely that 
the GP evolutionary stage is very brief. This 
is also to be expected from the enormous ener- 
getic release from the starburst that will quickly 
heat up and disperse the cold gas reservoir of 
GPs, quenching star-forming activities after a few 
Myr. The starburst nature of GP is also supported 
by STARLIGHT fits which suggest longer quiescent 
phases preceeding the current burst ranging from 
a few 10* yr in the case of GP232539 to -1 -10 
Gyr for GP004054 and GPl 13303. Note that the 
smoother SFH of GP232539, resulting into a sub- 
stantial intermediate-age stellar component is in- 
dicated by both RUN 1 and RUN 2. 

The overall evidence (stellar mass, SFH pat- 
tern and burst parameter) is consistent with the 
hypothesis of GPs being identifiable with previ- 
ous starburst-dominated phases in the lifetime of 
luminous BCDs in the local Universe. As appar- 
ent from HST/WFPC2 J?-band (F606W) images 
(Fig. 1), the GPs under study are as well com- 
patible to local BCDs with respect to their com- 
pactness (~5 kpc) and the irregular morphology of 
their high-surface brightness component. Further- 
more, similar to what is found for the main popu- 
lation of chemically evolved (7.6<12-F(0/H)<8.3) 
nearby BCDs (e.g., Papaderos et al. 1996; Cairos 
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et al. 2001, 2003; Nocskc ct al. 2003; Gil dc Paz & 
Madore 2005; Vaduvescu et al. 2006; Corbin et al. 
2006; Amorm et al. 2007, 2009), the GPs studied 
here show a more extended lower-surface bright- 
ness (LSB) envelope that is presumably due to an 
old stellar host. It is worth pointing out, on the 
other hand, that, contrary to most BCDs, the LSB 
host of GPs displays significant departures from 
ellipticity, suggesting a lesser degree of dynamical 
relaxation, or merger origin, which appears rather 
typical for the high-luminosity end of local BCDs 
(e.g. Telles et al. 1997; Ostlin et al. 2001; Bergvall 
& Osthn 2002; Adamo et al. 2011) This is particu- 
larly true for GP113303 for which HST data reveal 
significant LSB emission out to ~10 kpc southwest 
of its miclear region. 

The properties of the LSB component can be 
better quantified from the surface brightness pro- 
files (SBPs) in Fig. 8. The SBPs were derived 
using method iv in Papaderos et al. (2002) and 
were converted to the Vega system, following the 
prescriptions by Holtzman et al. (1995). They 
have been corrected for cosmological dimming as 
Mcorrcctod = ^observed " 10 log(l-l-z) and for Galac- 
tic extinction but no k corrections were applied. 
It can be seen that all SBPs show an outer expo- 
nential intensity decrease in their LSB envelope, 
as is typically the case for nearby BCDs. On the 
assumption that the LSB emission is of stellar ori- 
gin, one can infer from linear fits to the exponen- 
tial part of the SBPs the scale length and cen- 
tral surface brightness of CP hosts to 1.5-3 kpc 
and 19-21.5 mag arcsec"^, respectively. In Fig. 9 
we compare the structural properties of the host 
galaxy of GPs with those of local BCDs, dwarf ir- 
regulars (dis), LSB galaxies and dwarf ellipticals 
(dEs). It can be seen that GPs fall in the param- 
eter space that is populated by luminous BCDs, 
indicating that these two galaxy classes are very 
similar structurally. 

One word of caution is in order here. A roughly 
exponential intensity drop off is a generic prop- 
erty of extended nebular halos that are expected 
to arise in galaxies with strong starburst activity 
(Papaderos et al. 2002). An exponential profile 
in the LSB periphery of GPs should not therefore 
be taken as a foolproof sign of an extended stellar 
host galaxy. If the LSB envelope is partly or en- 
tirely due to extended nebiilar emission, then the 
absolute magnitude of the GP host will decrease 



by 0.75-1 mag, with a simultaneous brightening 
of the central surface brightness and decrease of 
the exponential scale length. This would move 
the GPs under study downstream in either plot, 
in better agreement to lower-luminosity BCDs. 

4.3. The nitrogen over-abundcince in GPs 

The good quality of the GTC data allowed us to 
derive the physical properties and chemical abun- 
dances of the ionized gas in the sample GPs with 
much better accuracy than previously done for 
these galaxies. 

The derived electron temperatures and densi- 
ties (Te([Oiii]) ~ 13000-15000 K, ne([Sii]) < 600 
cm~^) are consistent to those measured by AlO 
and by Izll for a larger sample of GPs. Moreover, 
these values are in general agreement with those 
typically measured in other samples of nearby 
star-forming dwarf galaxies (e.g., Campbell et al. 
1986; Masegosa et al. 1994; Lee et al. 2004; Hagele 
et al. 2008). 

The oxygen abmidances (12+log(0/H)) found 
for the sample GPs arc 7.81±0.14, 7.91±0.10, 
and 7.98±0.06 for GP232539, GP113303, and 
GP004054 respectively. These metallicity values 
are also consistent with those previously published 
using SDSS data (sec Table 1), and confirm that 
the three GPs are metal-poor galaxies. 

On the other hand, the three GP galaxies ap- 
pear to be relatively rich in nitrogen. The log- 
arithmic N/0 ratios determined for the sample 
GPs are -1.06 ± 0.04, -1.04 ± 0.08, and -0.99 
± 0.11 for GP004054, GP113303, and GP232539, 
respectively. These values are slightly higher (0.2 
dex on average) than those previously obtained 
by AlO. Differences can be attributed to large dif- 
ferences in S/N between GTC and SDSS spectra. 
Moreover, it is worth noting that in the calcula- 
tion of the nitrogen abundance we took into ac- 
count the [Nii]A6548A, which is clearly detected 
in the GTC spectra and mostly undetected in the 
SDSS spectra. The three GP galaxies studied 
here show N/0 ratios that are clearly higher than 
the usual upper limits of the plateau on the N/0 
vs. 0/H diagram and more compatible with those 
of emission-line galaxies with higher metallicities 
and similar sSFR (e.g., Mallery et al. 2007) for 
which secondary nitrogen production takes over. 
Using a large sample of galaxies that include the 
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GPs of C09, Izll noticed that for galaxies with 
EW(H;3)<100A there is a dependence of the N/0 
ratio on the total galaxy mass, implying some sec- 
ondary nitrogen production in the highest-mass 
galaxies (>1O^M0, see their Fig. 12). 

In order to explain the observed high N/0 ra- 
tios in metal-poor SFGs several mechanisms, in- 
cluding time delays on chemical enrichment, ef- 
fects of gas flows, and variations in the SFHs 
(e.g., recurrent bursts) have been proposed in the 
literature (see e.g., Garnett 1990; Pilyugin 1992; 
Koppen & Edmunds 1999; van Zee & Haynes 2006; 
MoUa et al. 2006; Pilyugin et al. 2012). Using a 
large sample of metal-poor emission-line galaxies, 
Izotov et al. (2006) found an apparent increase of 
N/0 with decreasing EW(H/3), best seen at inter- 
mediate metallicities. They interpreted this trend 
as an evidence of gradual enrichment of the Hii re- 
gions in nitrogen by massive stars from the most 
recent starburst. In this line, the discovery of WR 
stars in the three GPs studied here offer additional 
ingredients to the discussion. 

Localized chemical pollution by WR stars has 
been argued to explain the enhancement in nitro- 
gen abundance in the central starburst of nearby 
WR galaxies (e.g., NGC5253 Kobulnicky & Skill- 
man 1997; Lopez-Sanchez et al. 2007, 2011; James 
et al. 2009; Monreal-Ibero et al. 2010). These re- 
gions, with sizes of about tenths of parsecs usually 
show N/0 enhancements consistent with few WRs 
stars (Monreal-Ibero et al. 2010). In their inte- 
grated spectra, most WR galaxies with EW(H/3) 
< lOOA in Brinchmann et al. (2008) (see also 
Lopez-Sanchez & Esteban 2010) show higher N/0 
than non-WR galaxies, suggesting a rapid enrich- 
ment by WR stars. However, most of these WR 
galaxies with N/0 values comparable to our GPs 
show much higher oxygen, neon, and argon abun- 
dances at the same EW(H/3) (see their Fig. 15). 
Helium abundances in the three GP galaxies are 
roughly the same and not espcxdally high. Indeed, 
these values are comparable to those found in most 
BCDs (e.g., Pagel et al. 1992). As pointed out 
by Pagel et al. (1986, 1992) (see also Kobulnicky 
& Skillman (1997)) high He/H ratios are to be ex- 
pected in regions with strong WR pollution. How- 
ever, as remarked by Brinchmann et al. (2008), 
following results found in Galactic ring nebulae 
by Esteban et al. (1992), the expected average in- 
crease in He/H could be up to 0.7 dex lower than 



the expected increase in N/0. So the apparently 
normal He/H seen in the GPs may not give addi- 
tional clues in this respect. 

Finally, recent results from integral field spec- 
troscopy on three (HS 0128-h2832,HS 0837-1-4717 
and Mrk930) nitrogen overabundant BCDs by 
Perez-Montero et al. (2011), have shown that 
these galaxies display high N/0, which are con- 
stant at spatial scales of the order of several kpc. 
Interestingly, their N/0 values are rather similar 
to those of our target GPs. Perez-Montero et al. 
show the evolution of N/0 with time for spherical 
gaseous distributions of constant density with ini- 
tial 12-hlog(0/H)=8 and log(N/0)=-1.6 , which 
is polluted by the stellar winds coming from stel- 
lar clusters with Z = 0.004 and different stellar 
masses (from 10^ to lO^^M©). Importantly, their 
results show clearly that the number of WR stars 
estimated from spectroscopy cannot be responsi- 
ble for the observed enhancement in the N abun- 
dance across the observed spatial scales. Accord- 
ing to their results, in the case of the GPs it would 
be neccesary ionizing clusters with masses 2 or 3 
orders of magnitude higher to reach such degree 
of pollution in tlie spatial scale covered by our 
spectral aperture (see §3.3.1). This suggests that 
another chemical evolution scenario is required to 
explain the observed properties of the GPs. 

One possiblity is to think our galaxies as (chem- 
ically) evolved enough to have reach the secondary 
nitrogen regime. Koppen & Hensler (2005) used 
models to investigate the chemical evolution of 
galaxies during an episode of massive and sudden 
accretion of metal-poor gas due to, e.g., interac- 
tion with gas clouds or small gas-rich companions. 
In galaxies that has reached the secondary nitro- 
gen regime a rapid decrease of the oxygen abun- 
dance can be produced during the infall, remain- 
ing the N/0 unchanged, before a slower evolu- 
tion which leads back to the closed-box relation. 
This theoretical model and similar results from nu- 
merical simulations (e.g., Finlator & Dave 2008) 
led to AlO to suggest a scenario where the GPs 
would be experiencing the immediate effects of a 
recent, massive accretion of metal-poor gas (from 
the galaxy outskirt or beyond), possibly modu- 
lated by enriched outflows. This, at the same 
time, appears consistent with the SFH and ob- 
served morphology of the GP galaxies. 



13 



5. Summary and final remarks 

In this paper, wc have presented for the first 
time deep optical broad-band imaging and spec- 
troscopy for a sample of three GPs with a 10m- 
class telescope, the GTC. 

Chemical abundances and physical properties 
of the ionized gas were derived with high precision 
in the three galaxies, confirming previous findings. 
In particular, our results support several of the 
conclusions drawn by AlO: We find that these 
galaxies are metal poor (7.8<12-hlog(O/H)<8.0) 
just like nearby BCDs. However, in contrast to 
most BCDs, the studied GPs show a remarkably 
high N/0 ratio for their low oxygen abundance. 
This is in agreement with previous work by AlO 
(see also Pilyugin et al. 2012) who claimed after 
a re-analysis of SDSS data for the C09 sample of 
GPs that a significant fraction of these systems 
shows enhanced nitrogen abundances. 

The high S/N GTC spectra allowed to discover 
clear signatures of WR stars in the three galax- 
ies. One of the implications of this finding is that 
these WR stars might produce a localized chemi- 
cal pollution, and therefore explain the apparent 
nitrogen over- abundance of the galaxies. However, 
by comparing with similar results for BCDs pre- 
sented by Pcrcz-Montcro et al. (2011), the derived 
ratio between WR and O stars for GPs seem to 
be not enough to explain the apparent nitrogen 
over- abundance as due to local pollution by WRs. 

With the aim of gaining further insight into 
the SFH of GPs, we carried out a pilot spec- 
tral synthesis study of these systems using the 
code STARLIGHT in conjunction with the new- 
generation SSP library POPSTAR which includes 
nebular emission (runI). As a consistency check, 
we additionally employed STARLIGHT with purely 
stellar SSPs (run2), as well as an evolutionary 
synthesis code that self-consistently accounts for 
Balmer emission line equivalent widths (EWs) for 
a simplified SFH comprising an old and a young 
stellar population (run3). The central outcome 
from the STARLIGHT-I-POPSTAR model, qualita- 
tively confirmed through the RUN2&RUN3, is that 
the GPs under study currently undergo a major 
starbust producing between ^4% and 20% of their 
stellar mass. This result is supported by the de- 
tection, for the first time in these galaxies, of clear 
signatures of a significant number of WR stars in 



their spectra. The derived SFH of these systems 
shows large discontinuities between a few 10® yr 
and several Gyr, strongly suggesting that their for- 
mation history is dominated by starbursts. How- 
ever, as models imply, GPs are old galaxies having 
had formed most of their stellar mass several Gyr 
ago. The presence of an old stellar component 
could be directly spectroscopically verified in one 
of our GPs, adding further support to this con- 
clusion. Stellar masses derived from runI range 
between 1.7 x 10^ Mq and 2.7 x 10^ Mq, they 
are thus fairly comparable to those of luminous 
BCD galaxies in the nearby Universe. Moreover, 
a surface photometry study based on HST data in- 
dicates that our sample GPs possess an extended 
low-surface brightness (LSB) envelope that, if due 
to stellar emission, has structural properties com- 
patible to those of luminous BCDs in the nearby 
Universe. 

Clearly, to further elucidate the evolutionary 
scenarios of GPs we need a better understanding 
of their chemical enrichment and star formation 
histories. In particular, it would be important to 
compare accurate measurements of their chemical 
abundances, especially nitrogen, for a larger sam- 
ple of GPs with chemical evolution models tak- 
ing into account gas inflows and outflows. High- 
quality spectra with 8-lOm class telescopes have 
been proved to be indispensable to study in detail 
faint spectral features and to reduce uncertainties 
in the analysis and interpretation of the data. 
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Main properties of the galaxies from the literati.tre 
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Note. — Redshifts (z) were taken from SDSS DR7. E{B — V) values are from Schlegel et al. (1998). Magnitudes are not corrected from 
galactic and internal extinction. Absolute magnitudes were calculated as Mg = mg+ 5 — 5 log{Dj^)+ 2.5 log{l + z), where Dl is the 
luminosity distance and nig is the Sloan p'— band apparent magnitude {rripetro)- Half-light radius (re), obtained from WFPC2/HST optical 
broad-band imaging, and dust-corrected total SFRs are from Overzier et al. (2009). Logarithm of the FUV luminosity {Lpuy^ was taken 
from Hoopes et al. (2007) 

*From Overzier et al. (2009) 
''From Izotov et al. (2011) 

"^Values obtained and used by Amorm et al. (2010). they were derived using the direct method in the case of GP113303, and the N2 
calibration for GP232539. 



Table 2 

OSIRIS-GTC LOG OF observations 



ID 


CCD Chip 


Slit* /Filter 




Airmass 


Seeing 


Date 


Comment 


GP004054 


2 


Sloan z 


2250s 


1.10-1.25 


<1.0 


2010-07-08 


clear/dark 




2 


RIOOOB 


3870s 


1.1-1.2 


1.2 


2010-07-11 


spectroscopic-dust / dark 




2 


RIOOOR 


3870s 


1.13-1.26 


1.0 


2010-07-17 


clear/dark 




2 


RIOOOR 


3870s 


1.12-1.25 


1.0-1.1 


2010-08-02 


clear-dust / dark 


GP113303 


2 


Sloan z' 


2250s 


1.45-1.60 


1.0 


2010-05-21 


photometric / gray 




1 


RIOOOB 


3870s 


1.35-1.45 


<1.2 


2010-06-03 


photometric / gray 




1 


RIOOOR 


3870s 


1.80-2.00 


1.1-1.4 


2010-07-09 


clouds-dust / dark 


GP232539 


2 


Sloan z' 


2250s 


1.3-1.5 


<1.2 


2010-07-07 


clear-dust / dark 




1 


RIOOOB 


2580s 


1.14 


1.0 


2010-07-09 


cloudy-dust/ dark 




2 


RIOOOB 


1290s 


1.43 


1.4 


2010-07-12 


cloudy-dust /dark 




1 


RIOOOR 


3870s 


1.13-1.18 


1.4 


2010-07-10 


spectroscopic-dust / dark 



^In all cases the spectroscopic observations were done with a 0.8 arcsec wide slit located in the parallactic angle. 
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Table 3 

Measured and extinction-corrected line fluxes and line identifications 



A, Line id. 


f(A) 


GP004054 


GP 113303 


GP232539 






F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


3188 Hel 


0.47 


20 i 2 


28 it 4 






32 it 14 


51 i 22 


3499 Hel 


0.37 


1 1 4- 
±1 It z 


1 /I -I- 










3531 Hel 


0.36 




1 


25 it 2 


30 it 5 






3697 H17 


0.33 


Q 4- 1 


1 -1- 1 










3704 H16+Hel 


0.33 


I 1 4- 

I I It 


14 It 


10 It 


-Xl -X- A 
1 ( It 4 


fi -I- 1 
It 1 


I 1 -U 

II ± z 


3712 H15 


0.33 


12 i 1 


15 i 1 






10 it 5 


lO It i 


3727 [Oil] 


0.32 


1314 i 4 


iOoU it f 


A AO. -L 1/1 
144o it 14 


1 ft'Tc; -L 010 
iO/O it ALA 


1815 it 8 


A TO -L 111 
Z4 / y it 111 


3750 H12 


0.32 


23 i 2 


30 i 3 






9 it 2 


12 it 3 


3770 Hll 


0.31 


31 i 4 




zl i 


OA -4- 


14 i 


19 it 4 


3798 HIO 


0.31 


34 i 1 


43 i 3 


29 it 1 


33 it 4 


32 it 6 


43 it 8 


3820 Hel 


0.30 


i<5 It £1 


.^0 It 






ft -I- 1 
fS It 1 


1 1 -1- 


3835 H9 


0.30 


4- 
00 It £1 


( 1 It 4 


OJi It 


ftn -1- Q 

OU It 


OQ -1- Q 
AM It 


oy It 


3868 [Nelll] 


0.29 


QQQ 4- 1 n 
oOy It iU 


/ISft -1- Oft 
40D It ^0 


A '7C\ -\- 

4/U It 


EOQ _1_ ftO 
DOO It 


Qflfi -1- 1 
OUcS It LA 


/I no -1- OQ 

4uy It Ao 


3889 Hel+H8 


0.29 


1 KC\ 4- A 
iOU It 4 


1 Qft -1- 1 n 
loD It lU 


1 KA -\- A 
104 ± 4 


1 Tft -1- on 

1 / It A\J 


■\ AO -\- K 
140 It 


1 QQ -Lin 

icy It lu 


3916 Nil 


0.28 




it 


68 it 11 


78 it 15 






3968 [NoIII]+H7 


0.27 


287 It 4 


351 it 16 


252 it 6 


285 it 31 


182 it 6 


235 it 12 


4026 [NII]+HoI 


0.25 


OQ 4- /I 
ZO It 4 


28 it 5 






10 it 2 


13 it 2 


4102 H<5 


0.23 


247 i 3 


294 it 11 


2do it 10 


zyo it zy 


218 it 6 


272 it 11 


4243 [Foil] 


0.19 




1 

it 


25 it 3 


27 it 4 






4340 H7 


0.16 


463 It 4 


522 it 14 


418 it 4 


A A(\ -L OQ 


437 it 5 


508 it 13 


4363 [Olll] 


0.15 


4- Q 

/y It 


QQ -I- A 

oO It 4 


HA -\- Q 
<54 It 


on -1- ft 
yu It 


ft/1 -1-11 

D4 It 11 


'7 A -U 1 
f 4 It LA 


4471 Hel 


0.12 


Af\ 4- 1 
4:U It 1 


A A -1-1 
44 ± 1 


0/ It 


oy It 1 


00 It 


OT -X- A 
f It 4 


4571 Mgl] 


0.08 


( It i 


7-1-1 
/ It 1 










4658 [Folll] 


0.06 


18 i 2 


19 it 2 


18 it 4 


18 it 4 


18 it 2 


19 it 2 


4686 Holl 


0.05 


1 /I 4- 
±4t It z 












4702 [Felll] 


0.05 


3 ± 


it i 




1 

— it 


1 

— it 




4713 [ArIV]+HcI 


0.04 


12 i 1 


13 it 1 


-1- 
y It 


Q -1- 

y It ^ 


-I- 1 
It 1 


-I- 1 
It 1 


4740 [ArlV] 


0.03 


n 4- 1 

y it 1 


y it i 










4815 [Fell] 


0.01 


D it 1 


D it i 


16 it 5 


16 it 5 






4861 H/3 


0.00 


1 ono -1- 3 


1000 + '\ 


1000 + 6 


1 000 + 6 


1000 + 5 

—1— »j 


1000 + 5 


4881 [Felll] 


0.00 




- ± - 


19 it 2 


19 it 2 


7 it 1 


7 it 1 


4921 Hol 


-0.02 


13 ± 3 


13 ± 3 






4 it 2 


4 it 2 


4959 [Olll] 


-0.03 


2017 ± 8 


1978 ± 11 


1918 it 11 


1895 it 22 


1404 it 8 


1369 it 9 


4986 [Folll] 


-0.03 


14 ± 2 


13 ± 2 


14 ± 6 


13 ± 6 


12 ± 1 


12 ± 1 


5007 [OIII] 


-0.04 


6080 ± 13 


5907 ± 38 


5499 ± 27 


5404 ± 84 


4209 ± 24 


4058 ± 31 


5199 [Nl] 


-0.08 


12 ± 2 


11 ± 1 






30 ± 3 


28 ± 3 


5876 Hel 


-0.20 


131 ± 3 


112 ± 5 


128 it 2 


117 it 9 


141 it 5 


116 it 6 


6046 [Ol] 


-0.23 


31 ± 18 


26 ± 16 










6300 [OI] 


-0.26 


55 ± 2 


45 ± 3 


82 it 9 


72 it 11 


87 it 8 


68 it 7 


6312 [SIII] 


-0.26 


26 ± 6 


21 ± 5 


43 it 6 


38 it 7 


26 it 8 


20 it 6 


6364 [OI] 


-0.27 


15 ± 1 


12 ± 1 






28 it 7 


22 it 5 


6548 [Nil] 


-0.30 


64 ± 10 


51 it 8 


21 it 10 


18 it 9 


38 it 11 


28 it 9 


6563 Ha 


-0.30 


3588 ±12 


2861 it 138 


2952 it 33 


2573 it 301 


3878 it 49 


2907 it 125 


6584 [Nil] 


-0.30 


148 ± 6 


117 it 7 


172 ± 6 


150 it 20 


332 it 4 


248 it 11 


6678 Hoi 


-0.31 


41 ± 1 


32 ± 2 


141 it 10 


122 ± 17 


41 it 2 


30 it 2 


6717 [SII] 


-0.32 


193 ± 4 


152 it 8 


157 it 10 


136 it 19 


330 it 5 


243 it 11 


6731 [SII] 


-0.32 


149 ± 3 


117 it 6 


165 it 21 


142 it 25 


253 it 8 


185 it 10 


7065 Hel 


-0.36 


43 ± 8 


33 it 6 






37 it 13 


26 it 9 


7083 [Arl] 


-0.37 


8 ± 1 


6 it 1 






28 it 2 


20 it 2 


7136 [Arlll] 


-0.37 


91 ± 6 


69 it 6 






84 it 4 


59 it 4 


7388 [Fo II] 


-0.41 










17 it 1 


11 it 1 


7499 Hoi 


-0.42 


67 ± 2 


48 it 4 






13 it 2 


9 it 1 


c(H/3) 




0.33 ± 0.07 


0.20 ± 0.17 


0.42 ± 0.05 


-W(H/3) (A) 






156 




74 




60 


F(H/3) 






18.7 




15.4 




20.2 



Note. — Wavelengths in Col. 1 (A) arc rcst-framc. Measured [F(A)] and extinction-corrected [1(A)] emission-line 
fluxes are relative to F(H/3) = I(H/3) = 1000 with their corresponding errors. We also give for each galaxy the 
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corresponding constants of reddening [c(H/3)], the H/3 emission-line extinction-corrected total flux in units of 10 
erg cm^ A~^, and the H/9 equivalent width. 
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Table 4 

Derived physical conditions and chemical abundances 



Object 


GP004054 


GP 113303 


GP232539 


ne([Sll]) (cm-3) 


< 260 


< 680 


< 210 


te([Om]) (10^ K) 


1.34 ± 0.02 


1.39 ± 0.04 


1.46 ± 0.11 


te([Oll])^ (10* K) 


1.52 ± 0.02 


1.60 ± 0.02 


1.60 ± 0.07 


te([Nii]) (10-1 K) 


1.26 ± 0.01 


1.29 ± 0.02 


1.32 ± 0.04 


te([SlIl])b (10-* K) 


1.27 ± 0.15 


1.34 ± 0.16 


1.42 ± 0.20 


y+ (4471) 


0.090 ± 0.002 


0.078 ± 0.014 


0.077 ± 0.009 


y+ (5876) 


0.088 ± 0.004 


0.086 ± 0.007 


0.090 ± 0.004 


y+ (6678) 


0.088 ± 0.006 


- 


0.084 ± 0.007 


y+ (7065)'= 


0.110 ± 0.020 




0.082 ± 0.028 


y2+ (4686) 


0.0012 ± 0.0001 




- 


y 


0.089 ± 0.002 


- ± 0.004 


0.087 ± 0.004 


12+log(0+/H+) 


7.19 ± 0.04 


7.15 ± 0.07 


7.30 ± 0.10 


12+log(0+2/H+) 


7.91 ± 0.02 


7.82 ± 0.03 


7.65 ± 0.05 


12+log(0/H) 


7.98 ± 0.06 


7.91 ± 0.10 


7.81 ± 0.14 


12+log(N+/H+) 


6.13 ± 0.01 


6.11 ± 0.02 


6.31 ± 0.04 


ICF(N+) 


6.18 ± 0.18 


5.73 ± 0.23 


3.21 ± 0.59 


12+Iog(N/H) 


6.92 ± 0.02 


6.87 ± 0.03 


6.82 ± 0.08 


log(N/0) 


-1.06 ± 0.04 


-1.04 ± 0.08 


-0.99 ± 0.11 


12+log(Nc+VH+) 


7.27 ± 0.03 


7.26 ± 0.06 


7.08 ± 0.08 


ICF(Nc+^) 


1.08 ± 0.04 


1.08 ± 0.06 


1.11 ± 0.11 


12+Iog(Ne/H) 


7.30 ± 0.03 


7.30 ± 0.06 


7.13 ± 0.09 


log(Ne/0) 


-0.68 ± 0.07 


-0.61 ± 0.12 


-0.68 ± 0.17 


12+log(Ar+^/H+) 


5.54 ± 0.04 




5.40 ± 0.04 


12+log(Ar+3/H+) 


4.57 ± 0.05 






ICF(Ar) 


1.04 ± 0.04 




3.00 ± 1.44 


12+log(Ar/H) 


5.60 ± 0.10 




5.88 ± 0.22 


log(Ar/0) 


-2.31 ± 0.10 




-1.77 ± 0.23 


12+log(S+/H+) 


4.74 ± 0.15 


4.76 ± 0.25 


4.94 ± 0.14 


12+log(S+VH+) 


6.32 ± 0.25 


6.49 ± 0.24 


6.13 ± 0.27 


ICF(S) 


1.51 ± 0.40 


1.55 ± 0.55 


1.27 ± 0.27 


12+log(S/H) 


6.51 ± 0.35 


6.69 ± 0.62 


6.26 ± 0.36 


log(S/0) 


-1.48 ± 0.36 


-1.22 ± 0.43 


-1.55 ± 0.39 


12+log(Fe^+ 


5.69 ± 0.05 


5.62 ± 0.09 


5.59 ± 0.08 


ICF(Fc2+ 


6.3 ± 4.8 


7.0 ± 6.0 


3.7 ± 2.7 


12+log(Fe/H) 


6.49 ± 0.25 


6.39 ± 0.36 


6.11 ± 0.29 


log(Fe/0) 


-1.49 ± 0.26 


-1.52 ± 0.38 


-1.69 ± 0.32 



''Derived from te([Oin]) using the models described in Perez-Montero & 
Diaz (2003) and Perez-Montero & Contini (2009) for te([On]) and te([Nil]), 
respectively 

''Derived from te([Oni]) using the empirical relation described in Hagele 

et al. (2006) 

■^Hel 7065A should be affected by radiative transfer effects. However, 
we considered here their measurements in our calculations since ?/+(7065) 
values are, within the uncertanties, in agreement with the values derived 
from the other helium lines. 
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Fig. 1.— Green Pea imagery: GTC-OSIRIS z'-band images of GP004054 (upper), GP 113303 (middle), 
and GP 232539 (bottom), are shown in surface brightness grey scale, as indicated by the bar in the lower 
right panel. In all cases contours indicate the fiz =26 mag arcsec"^ level. HST WFPC2 F606W (R) images, 
spatially resolving the central regions of these galaxies, are showed as insets in their lower left corner. North 
is up and east to the left. 
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Fig. 2. — The Deep GTC-OSIRIS spectra for the three GPs under study. In order to highhght faint emission 
lines each spectrum is zoomed in the lower panel. 
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4400 4600 4800 5000 



Fig. 6. — Spectra of the galaxies zoomed in the 
region around the Wolf-Rayet blue bump. 




5000 5100 5200 5300 



Fig. 7. — Spectrum of GP113303 zoomed in the 
region around the Mgi absorption line after con- 
tinuum subtraction. The central wavelength of the 
Mgi A5173 is indicated by the vertical dashed line. 
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Table 5 

Model results from starlight (run 1 and run 2) 



Galaxy nebular continuum'^ (%) 
[O II] H/3 Ha 


, young 

(%) 


M,, total 
(log Mq) 




RL'K 1; POPSTAR MODELS 










GP004054 10 18 


35 


22 


9.18 


1.05 


GP113303 3.5 6 


11 


4.2 


9.44 


0.86 


GP232539 8 10.5 


20 


11.3 


9.38 


1.54 


RUN 2; B&C MODELS 










GP004054 




15.7 


9.22 


1.17 


GP113303 




1.9 


9.68 


0.95 


GP232539 




5.4 


9.90 


1.37 



^For three difi'ereni wavelength regions (around [On], HP, and Ha) , 
columns 2, 3, and 4 give the estimated fraction nebular continuum emission 
derived from starlight fits to each galaxy spectrum using popstar models 
(see text) 



Table 6 

Results from the two-component evolutionary synthesis model (run 



Galaxy 


^burst 


bpar 


C(H/3) 


EW(Hc«) 


EW(H/3) 


, young 






(Myr) 






(A) 


(A) 


(%) 




GP004054 


4 


5.3 


0.41 


894 


139 


8.3 


1.5 


GP113303 


5 


4.2 


0.32 


445 


76 


6.7 


2.1 


GP232539 


5 


2.8 


0.43 


344 


53 


4.4 


1.4 



Note. — Column 2: age of the young stellar component. Column 3: burst parameter. 
Column 4: reddening constant predicted by the evolutionary model. Columns 5 and 
6: predicted equivalent widths of Ha and H/3 emission lines. 



Table 7 

Wolf-Rayet features detected in the galaxy sample 



Galaxy 


/bb(A)'' 


logiBB*' 


EWbb 


N(WR,)'= 




GP004054 


34±4 


40.1 


4.2±0.5 


806±97 


O.OllibO.OOl 


GP113303 


8a±15 


40.8 


5.5±1.0 


1239±228 


0.029±0.005 


GP232539 


36±6 


40.3 


1.9±0.3 


1145±191 


0.012±0.002 



^Extinction-corrected, emission-line flux measured for the WR blue bump 
relative to I(H/3) — 1000 with their corresponding errors. 

^uminosity of the WR blue bump calculated using /bb(A) and the dis- 
tances adopted throughout the paper 

"^Number of WR stars as derived from Lbb and adopting a mean luminosity 
for a WR star logLwR* = 37.248 from starburst99 models with stellar 
metallicity Z = 0.004 (Perez-Montero et al. 2010) 

*^WR to O star number ratio and its errors. The number of O stars 
were derived from Lhck and adopting a mean Ha luminosity for an O star 
logLHa,o* = 37.21. 
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Jl 13303 



2 4 6 8 

R (kpc) 

Fig. 8. — Surface brightness profiles (SBPs) of 
the three GPs derived from archival HST/WFPC2 
R— band (F606W) images. The SBPs have been 
corrected for Galactic extinction and cosmological 
dimming but no k corrections were applied. The 
SBP of J232539 has been shifted by -1 mag for the 
sake of clarity. It can be seen that all SBPs show 
an extended exponential LSB envelope dominat- 
ing for /i > 22 — 23 mag arcsec"^. Linear fits to 
the LSB component (gray lines) yield scale lengths 
and central surface brightness levels of 1.5-3 kpc 
and 19-21.5 mag arcsec"^, respectively. 




Fig. 9. — Comparison of the structural proper- 
ties of the LSB component of our sample GPs 
(large crosses) with those of other types of low- 
mass galaxies compiled in Papaderos et al. (2008). 
Photometric quantities derived from fitting of R 
band profiles in Fig. 8 have been converted into 
B band assuming a B — R — 1 mag. The lower 
and upper panel compare, respectively, the central 
surface brightness and the exponential scale length 
vs the absolute magnitude of the LSB component 
of GPs with those of BCDs, dwarf irregulars (dis), 
dwarf ellipticals (dEs) and LSB galaxies. It can be 
seen that GPs occupy roughly the same parameter 
space as luminous BCDs. 
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wavelength (A) 



log (t [yr]) 



Fig. 3. — upper panel: Best-fitting synthetic SED based on popstar SSPs (runI; blue color), overimposed 
on the rest-frame observed spectrum of GP004054 (orange), normalized at 4170 A. The smaller plots on the 
right show the luminosity and mass contribution of individual SSPs (upper and lower panel, respectively). 
The age distribution of the library SSPs is illustrated by thin vertical lines in the upper panel, lower panel: 
Fit to the observed spectrum based on purely stellar SSPs from Bruzual & Chariot (run2). The color coding 
is as in the upper large panel. Vertical strips mark regions that have been flagged prior to spectral fitting. 
The best-fitting synthetic SEDs from a two-component evolutionary synthesis model (run3) that comprises 
an old and a young stellar component (red and blue shaded area, respectively) , and self-consistently accounts 
for the observed Balmer Ha and H/3 EWs, is overlaid in black color. 
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wavelength (A) log (t [yr]) 

Fig. 4.— Same as Fig.3 for GP113303. 
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